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摘 要 :本 研究 建立 了 纤维 增强 复合 材料 层 合板 混合 模式 弯曲 破坏 试 件 有 限 元 模型 ,采用 内 聚 力 模 
型 和 黏 结 单元 模拟 复合 材料 界面 层 的 力学 行为 ,对 现 有 复合 材料 层 合 板 混合 模式 分 层 破 坏 试验 进 
行 数值 模拟 ,预测 与 试验 荷载 -位 移 曲 线 吻 合 良好 。 根 据 已 建立 的 有 限 元 模型 ,探讨 了 内 聚 力 模型 
的 材料 参数 如 界面 单元 刚度 .界面 破坏 强度 .断裂 能 以 及 单元 尺寸 .不 同 材料 软化 法 则 和 材料 破坏 
准则 等 对 弯曲 试 件 混合 模式 破坏 试验 数值 模拟 准确 性 和 收敛 性 的 影响 。 结 果 表 明 ,合理 选择 单元 
尺寸 可 保证 数值 模拟 的 精确 性 ;合理 选用 材料 软化 法 则 和 界面 单元 刚度 对 运算 效率 提高 尤为 关键 ; 
断裂 能 取 值 必须 足够 准确 以 获得 合理 的 数值 分 析 结 果 ; 界 面 极限 强度 的 选择 对 运算 结果 影响 微小 。 
而 材料 破坏 准则 的 选取 对 混合 模式 分 层 破坏 数值 模拟 运算 效率 区 别 明 显 。 研 究 结果 为 更 有 效 地 使 
用 内 聚 力 模型 及 黏 结 单元 提供 了 参考 指导 作用 。 有 利于 工程 师 在 进行 复合 材料 层 合板 分 层 失 效 分 
析 遇 到 收敛 困难 时 ,更 好 地 选择 相关 参数 和 准则 以 保证 数值 计算 的 收敛 性 ,精确 性 和 运算 效率 。 
关键 词 : 复 合 材料 ; 内 聚 力 模型 ; 黏 结 单元 ; 分 层 破坏 

中 图 分 类 号 :TB332 ”文献 标志 码 :A DOI:10. 11776/j. issn. 1000-4939. 2023. 03. 006 


Numerical simulations of mixed-mode delamination 
in fiber-reinforced composite materials 
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Mechanics and Construction Engineering, Jinan University ,510632 Guangzhou , China) 


Abstract; A finite element ( FE) model is developed for simulating mixed-mode-bending delamination in 
fiber-reinforced composite laminates. A built-in cohesive zone model based on cohesive elements of the fi- 
nite element procedure ABAQUS V6. 14 is adopted to simulate the mechanical behavior of interfaces be- 
tween adjacent composite layers. Numerical simulations of mixed-mode-bending specimens were performed 
using the developed FE model. The predicted load versus displacement curves agree well with their experi- 
mental counterparts. The effects of the cohesive zone model parameters and criteria , such as the stiffness of 


cohesive elements , finite element size , interface failure strengths , fracture energies , and mixed-mode delam- 
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ination softening law and propagation criterion on the accuracy and convergence of numerical simulations 


of the mixed-mode-bending specimens , are discussed. The following conclusions can be drawn: proper se- 


lection of finite element size can ensure the accuracy of numerical simulations; appropriate choices of the 


interface cohesive element stiffness and mixed-mode delamination softening law are critical to improving 


the efficiency of numerical simulations; fracture energy should be sufficiently accurate for ensuring good 


numerical results; the choice of the interface failure strength has little effects on the simulation results. The 


choice of the delamination propagation criterion has distinct effects on the efficiency of numerical simula- 


tions of the mixed-mode delamination. The research provides a good reference for more effective use of co- 


hesive zone model and cohesive elements. It is helpful for engineers to select relevant cohesive zone model 


parameters and criteria to ensure the accuracy , computational convergence and efficiency of numerical a- 


nalysis when encountering computational difficulties in delamination analysis of composite materials. 
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纤维 增强 复合 材料 具有 高 比 模 量 、 比 强度 以 及 
秆 画 刚度 和 强度 可 设计 等 优点 ,被 广泛 应 用 于 土木 
下 程 .机械 工程 .海洋 工程 和 军事 工程 等 工程 领域 ， 
机 有 良好 的 发 展 前 景 。 但 是 与 纤维 相 比 ,基体 的 强 


许 相 对 较 弱 , 使 得 复合 材料 的 层 间 破 坏 模式 即 分 层 
看 环 成 为 了 最 常见 .最 主要 的 破坏 形式 。 复 合 材 
立 蒋 分 层 破坏 一 直 是 复合 材料 领域 的 研究 重点 。 

之 在 实际 工程 中 ,复合 材料 承受 复杂 的 荷载 ,混合 
模 瓶 分 层 破坏 更 为 普遍 。 国 内 外 学 者 对 复合 材料 混 
从 模式 分 层 破坏 进行 了 大 量 的 研究 , BENZEGGAGH 
等 (采用 试验 方法 研究 了 复合 材料 混合 模式 下 的 分 
Vc I E , 研究 表明 混合 模式 弯曲 破坏 试 件 
(fnixed-mode bending apparatus, Jr; X fj fy MMB ) 可 
以 很 好 地 模拟 混合 模式 分 层 破坏 ,而 且 数值 模拟 的 
荷 虐 -位 移 曲线 的 形状 与 混合 比例 (mixed-mode rati- 
oXG,/ G, (其 中 , G1 为 了 型 单一 模式 破坏 时 的 能 量 
释放 率 ，G 为 模型 在 当前 荷载 下 分 别 在 型 、 开 型 
和 亚 型 分 层 破坏 中 的 能 量 释放 率 之 和 ) 的 取 值 有 关 。 
MIYAGAWA 等 中 利用 拉 曼 光谱 法 测定 L, 
层 破坏 层 间 裂纹 尖端 附近 的 应 变 分 布 ,计算 得 出 的 
断裂 彻 性 与 测定 值 基本 一 致 。 庄 蔚 敏 等 通过 三 点 
弯曲 试验 ,研究 碳纤维 复合 材料 在 弯曲 载荷 作用 下 
的 失效 形式 与 损伤 过 程 ,发 现 碳纤维 复合 材料 层 合 


生 、 损 伤 的 积累 以 及 材料 的 最 终 失 效 破坏 做 了 综合 
性 的 回顾 。WHITNEY 等 以 及 KIM 等 中 采用 强度 
理论 预测 模型 来 判断 界面 是 否 发 生 分 层 , 当 材 料 内 
部 某 点 处 的 应 力 或 一 定 区 域 的 平均 应 力 超过 界面 强 
度 时 即 判 断 界面 发 生 分 层 损伤 。ALLIX 257 
CORIGLIANO 等 (中 针对 层 合板 界面 层 建立 了 损伤 
模型 来 预测 界面 的 分 层 损 伤 。 内 聚 力 理论 最 早 由 
DUGDALE00 和 BARENBLATT " 1 提出 ,用 于 描述 材 
料 界 面 分 离 时 的 断裂 过 程 。DUCDALE'"1 对 含 穿 透 
裂纹 的 薄 钢 板 进行 拉 伸 试验 ,发 现 裂 纹 尖 端 塑性 变 
形 区 为 扁平 带 状 区 域 , 并 将 其 简化 为 内 聚 区 。BA- 
RENBLATT'2 在 对 脆性 材料 断裂 过 程 的 研究 中 , 认 
为 内 聚 区 的 应 力 分 布 是 关于 该 点 与 裂纹 尖端 之 间距 
离 的 函数 。 为 准确 计算 该 区 域内 的 应 力 ,大 部 分 研 
究 引 入 内 聚 区 长 度 的 概念 '。 内 聚 区 长 度 的 定义 
HH ,从 完全 破坏 失去 承载 能 力 的 黏 结 单元 到 应 力 
达到 最 大 值 的 黏 结 单元 之 间 的 距离 。 在 内 聚 区 中 ， 
应 力 梯 度 变化 比 其 他 部 位 大 得 多 。 随 着 计算 机 技术 
的 发 展 , 许 多 学 者 采用 有 限 元 程序 材料 库 中 的 内 聚 
力 模 型 和 笑 结 单元 的 方法 “中研 究 分 层 损 伤 问题 。 

HUI 45^ 通过 试验 发 现 当 材料 弹性 模 量 小 于 
界面 茜 结 强度 时 ,界面 分 层 裂 纹 的 发 展 将 发 生 钝 化 
现象 ,并 通过 有 限 元 方法 采用 内 聚 力 单 元 对 裂纹 发 


板 弯曲 损伤 过 程 中 , 层 间 裂 纹 和 层 内 裂纹 扩展 速度 
不 同 , YE A Az He 2T Af IBI ER HH E Ap RES R, Ji pe] n 
纹 已 经 开始 出 现 。 云 新 部 等 ”对 斜纹 编织 碳纤维 / 
环 氧 树 脂 复合 材料 型 分 层 开裂 性 能 进行 了 测试 ， 
发 现 编织 纤维 增强 复合 材料 分 层 破 坏 行为 除 层 间 开 
袭 外 还 存在 纬 向 纤维 脱 务 。BOLOTIN' 和 根据 裂纹 在 
层 合板 的 位 置 不 同 , 对 层 内 及 近 表 面 分 层 裂 纹 的 产 


展 进 行 模拟 ,解析 解 与 数值 解 吻合 良好 。YANC 
等 ”建立 了 一 个 内 聚 力 模型 模拟 复合 材料 层 合板 
的 三 维 分 层 损伤 ,发 现 与 虚拟 裂纹 闭合 法 相 比 ,内 聚 
力 模型 能 正确 地 预测 试 件 的 分 层 及 辟 裂 裂纹 形状 ， 
但 是 内 聚 力 模型 的 相关 参数 需要 进一步 确定 。 
XLE 对 ASA/PEEK 复合 材料 1、 本 型 和 混 
合 模式 分 层 破 坏 进 行 有 限 元 数值 模拟 ,分 别 采 用 朝 
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第 3 期 
指数 准则 77 30 BK 准则 "作为 界面 失效 判 据 ,发 现 
采用 BK 准则 比 窜 指 数 准 则 预测 得 到 的 奏 载 -位 移 曲 
线 与 试验 曲线 更 为 接近 ,BK 准则 中 与 材料 相关 的 常 
数 的 选取 对 数值 模拟 从 载 峰值 的 精度 有 一 定 影 
响 ; 当 MMB 试验 破坏 模式 混合 比例 为 50% 时 ,7 取 
1 时 数值 模拟 荷载 峰值 与 试验 值 更 为 接近 。 国 内 学 
者 肖 先 等 分别 采用 BK 准则 和 短 指 数 准则 对 碳 纤 
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1 内 聚 力 理论 模型 


本 研究 采用 有 限 元 程序 ABAQUS V6. 14 材料 库 
中 的 内 聚 力 模型 和 COH3D8 黏 结 单元 模拟 复合 材料 
层 合 板 加 载 过 程 中 的 分 层 破坏 行为 ,可 同时 模拟 分 
层 破坏 的 发 生 与 发 展 过 程 。 通 常 ,在 复合 材料 单 层 


维 增强 树脂 基 材 料 的 单 向 层 合板 混合 模式 分 层 扩展 
进行 数值 模拟 , 发 现 采用 模型 参数 y = 1 n deis 
数 准 则 计算 得 到 的 混合 模式 断裂 能 释放 率 数值 仿 
大 ,而 采用 模型 参数 7 =2 时 BK 准则 计算 得 到 的 混 
合 模式 断裂 能 释放 率 与 理论 预测 值 更 吻合 ;SKOVA- 
JSA 等 中 针对 I 型 分 层 破坏 ,对 数值 模拟 与 解析 计 
算 进行 了 对 比 ,发 现 采用 网 格 尺寸 为 1 mm 或 2 mm 
的 线性 六 面体 单元 进行 数值 模拟 可 得 到 足够 精确 的 
结果 ; CAMANHO 等 5 发现, 宕 指数 准则 中 参数 y 
联 士 更 适合 用 于 模拟 ASA/PEEK 纤维 增强 复合 材料 
涯 全 模式 下 的 分 层 破坏 。 

一 采用 有 限 元 程序 中 的 内 聚 力 模型 和 黏 结 单元 模 
拟 复合 材料 的 分 层 破坏 力学 行为 是 普遍 的 做 法 。 复 
低 树 料 的 实际 结构 应 用 中 ,由 于 受 力 复杂 ,混合 模式 
和 处 刘 破 坏 的 发 生 更 为 普遍 。 因 此 ,对 复合 材料 的 分 
属 生 坏 模 拟 ,应 采用 考虑 混合 模式 破坏 的 内 聚 力 模 
型 进行 模拟 ,而 不 能 仅 采 用 单一 模式 的 内 聚 力 模型 。 
一 对 夭 结 单元 的 属性 定义 包含 本 构 模 型 ( 即 内 聚 
力 异 型 ) 以 及 性 能 参数 两 部 分 ,其 中 本 构 模型 包含 材 
料 矣 伤 发 生 准则 、 材 料 软 化 法 则 、 材 料 破坏 准则 ;性 
能 窖 数 包括 界面 刚度 .界面 强度 .断裂 能 。 以 上 准则 
及 容 数 共同 描述 黏 结 单元 的 力学 行为 ,要 使 用 夭 结 
单元 准确 地 模拟 复合 材料 分 层 破坏 过 程 ,必须 合理 
地 设置 上 述 各 个 准则 和 参数 。 

上 述 工作 大 部 分 是 单独 探讨 某 个 参数 对 分 层 破 
坏 数值 模拟 结果 的 影响 。 尚 未 发 现 有 研究 针对 基于 
不 同 材料 软化 法 则 、 材 料 破坏 准则 的 内 聚 力 模型 和 
不 同 模型 参数 如 界面 刚度 、 界 面 强度 等 对 复合 材料 
混合 模式 界面 分 层 破坏 进行 参数 化 分 析 , 探 讨 这 些 
因素 对 数值 模拟 结果 的 影响 ,因此 有 必要 开展 本 研 
究 工作 。 

本 研究 建立 纤维 增强 复合 材料 混合 模式 弯曲 破 
坏 试 件 有 限 元 模型 ,并 考虑 上 述 各 因素 对 数值 模拟 
结果 的 影响 ,可 为 复合 材料 结构 工程 设计 人 员 或 研 
究 人 员 更 有 效 地 利用 内 聚 力 模型 及 香 结 单元 模拟 复 
合 材 料 分 层 破坏 行为 提供 参考 指导 意见 。 


间 放 置 一 层 务 结 单元 模拟 界面 层 。 本 节 介 绍 内 聚 力 
模型 中 的 牵引 力 - 位 移 关 系 、 界 面 刚 度 、 材 料 损 伤 发 
生 准 则 、 材 料 软 化 法 则 、 材 料 破坏 准则 等 。 


1.1 牵引 力 -位 移 关 系 


黏 结 单元 的 本 构 关 系 由 牵引 力 -位 移 关 系 定义 。 
本 构 关 系 曲 线 分 为 2 段 , 第 一 段 通 常 采 用 线 弹 性 的 
牵引 力 -位 移 关 系 描述 分 层 损 伤 发 生前 界面 层 的 力 
学 响应 ,第 二 段 则 描述 分 层 损伤 发 展 后 的 材料 软化 
行为 。 当 分 层 损 伤 发 生 后 ,材料 属性 开始 退化 ,界面 
层 的 损伤 演化 由 材料 软化 法 则 控制 。 务 结 单元 的 牵 
引力 -位 移 关 系 可 表示 为 

ô, 


T, K, 0 0 
T, |= | 0 K 0 | 
T, 0 0 天 八 6， 


其 中 :zs 代 表 以 黏 结 单元 中 面 中 心 点 为 原点 的 局 
部 直角 坐标 系 的 3 个 方向 ,n 为 沿 单元 厚度 方向 ,对 
应 于 工 型 分 层 破坏 ,s 和 + 方向 代表 单元 中 面 平行 于 
纤维 和 垂直 于 纤维 的 方向 ,对 应 于 开 型 和 亚 型 剪 切 
分 层 破坏 ; o, 是 平面 外 牵引 力 ; r, 和 7, 分 别 为 在 垂 
直 于 单元 厚度 方向 的 平面 上 平行 于 纤维 方向 (界面 
内 第 一 剪 切 方向 ) 和 垂直 于 纤维 方向 (界面 内 第 二 
剪 切 方向 ) 的 剪 切 牵引 力 ; 6，、5，、 5 RI K, . K, K, 
为 相应 的 相对 位 移 和 刚度 。 刚 度 的 定义 如 下 


8, 
(1) 


K? , ô; x 
K =l -dK . (2) 
0 , 6,26 


其 中 :i=n,s,t; K? 是 初始 刚度 ; 0; 405; 分 别 是 每 种 
破坏 模式 下 材料 分 层 损伤 开始 发 生 以 及 材料 完全 破 
坏 时 对 应 的 相对 位 移 ; d; 为 描述 分 层 破坏 发 展 的 损 
伤 变量 ,其 计算 方式 与 选用 的 材料 软化 法 则 相关 。 


1.2 界面 刚度 


在 分 层 破 坏 发 生前 ,实际 上 界面 层 并 不 会 发 生 
变形 ,而 由 于 牵引 力 -位 移 关 系 中 引入 了 有 限 的 刚 
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度 , 导 致 分 层 破 坏 时 界面 层 会 发 生 拉 伸 变 形 ; 同 时 由 
于 梁 截 面 发 生 普 井 , 有 些 黏 结 单元 还 会 出 现 负 法 向 
应 力 和 压缩 变形 ,与 实际 情况 不 符合 。 因 此 ,界面 刚 
度 的 设置 应 尽量 保证 其 变形 与 实际 相符 。 当 设置 的 
刚度 较 小 时 ,条 结 单元 层 产生 较 大 的 相对 位 移 ,变形 
相对 于 结构 整体 变形 不 可 忽略 ,结构 受 力 及 变形 与 
实际 变形 发 生 偏离 ,将 得 到 错误 的 数值 模拟 结果 ; 当 
刚度 设置 过 大 时 , 则 可 能 出 现 计算 困难 而 使 数值 模 
拟 的 计算 效率 变 低 。 因 而 要 选择 足够 大 的 刚度 值 以 
保证 相 邻 复合 材料 单 层 具有 足够 坚固 的 连接 ,又 必 
须 在 合理 的 范围 内 从 而 避免 数值 计算 上 的 问题 , 提 
国内 外 学 者 提出 了 若干 种 确定 刚度 值 大 小 的 方 
法 。CAMANHO 等 “用 刚度 值 为 10 N/mm 模拟 石 
墨 = 环 氧 树脂 复合 材料 分 层 破 坏 并 得 到 了 准确 的 结 
SE-ZOU 等 25 建议 合适 的 刚度 值 范围 应 是 10 ~ 
1[_ 倍 的 界面 层 剪 切 或 者 抗 拉 强度 值 。TURON A 
"ACRI RE f c] 各 影响 复合 材料 整体 有 效 弹 性 


模 量 ， 
ee 2 (3) 
i Kt 


jp. nv 性 模 量 ; E, NES 
构 漳 层 厚度 方向 的 弹性 模 量 ;K 为 界面 层 刚度 ;t 为 
复 众 材料 层 合板 厚度 。 


ono, 当 P eK IM 


吉 构 的 有 效 弹 性 模 
E, 


ALETE $K-= , 当 


B Ho ,由 于 界面 层 — 失 小 
于 5% ,可 保证 界面 层 刚 度 对 结构 有 效 弹 性 模 量 的 
影响 足够 小 ,使 得 数值 模拟 足够 准确 。 在 本 研究 后 
续 的 参数 化 数值 模拟 中 ,其 中 采用 的 两 个 刚度 值 即 


di Kk = ÊF 分 别 取 且 = 20,50 计算 得 到 。 


1.3 材料 损伤 发 生 准 则 


为 预测 层 间 分 层 损伤 的 初始 ,采用 损伤 发 生 准 
则 进行 判断 。 当 黏 结 单元 的 应 力 状 态 满足 一 定 的 条 
件 时 ,分 层 破坏 即 开 始 发 生 , 同 时 材料 属性 开始 退 
化 。 复 合 材料 层 合板 分 层 破 坏 失 效 模 式 包括 工 型 、 
开 型 . 亚 型 和 混合 模式 分 层 失效 ,针对 不 同 模式 的 分 
层 损伤 起 始 , 需 要 用 不 同 的 材料 损伤 发 生 准 则 进行 
判断 。 在 ABAQUS V6. 14 有 限 元 程序 中 ,提供 了 4 
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种 损伤 发 生 准 则 :最 大 名 义 应 力 准 则 .最 大 名 义 应 
准则 、 二 次 名 义 应 力 准则 、 二 次 名 义 应 变 准则 。 
对 于 单一 模式 的 分 层 失 效 , 可 通过 比较 界面 应 


力 或 强度 来 判断 界面 是 否 发 生 分 层 损伤 ,一 般 采 用 
最 大 名 义 应 力 准 则 作为 判断 依据 ,其 表达 式 如 下 ， 
max [2 E 1 (4) 
Ca T7, 7, 


其 中 : o%、7，、7 分 别 代 表 黏 结 单元 变形 沿 界面 法 
线 方向 、 界 面 内 第 一 剪 切 方向 和 第 二 剪 切 方向 时 牵 
引力 -相对 位 移 曲 线 的 峰值 ,也 即 该 方向 务 结 单元 的 
极限 强度 ; x» = C |æ] +x)/2。 

对 于 混合 模式 分 层 失效 ,分 层 可 以 发 生 在 任何 
一 种 模式 的 层 间 应 力 达 到 它 的 极限 强度 之 前 。 
此 ,不 能 采用 最 大 名 义 应 力 准 则 或 最 大 名 义 应 变 准 
则 来 判断 混合 模式 的 分 层 失 效 。 

在 混合 模式 分 层 破 坏 中 ,一 般 采 用 二 次 名 义 应 
力 准则 或 二 次 名 义 应 变 准 则 来 判断 复合 材料 层 合板 
是 否 发 生 分 层 损伤 ,其 中 二 次 名 义 应 力 准则 表达 式 
如 下 。 


[22] «[5] [8] =e»p<1 (5) 
其 中 , 当 ew =1 时 ,表明 材料 发 生 了 分 层 损 伤 。 
而 二 次 名 义 应 变 准 则 表达 式 ， 


《En7 el &l 
| A ji + is] + is] e; «1 (6) 
其 中 : ea、e，、eo 分 别 代 表 黏 结 单元 变形 沿 界 面 法 
线 方向 .界面 内 第 一 剪 切 方向 和 第 二 剪 切 方向 时 黏 
结 单元 的 极限 应 变 ; es，、s，、e, 代表 当前 状态 界面 
法 线 方向 .界面 内 第 一 剪 切 方向 和 第 二 剪 切 方向 的 
应 变 ; 当 e, =1 时 ,表明 材料 发 生 了 分 层 损伤 。 本 研 


究 在 混合 模式 弯曲 破坏 数值 模拟 中 采用 了 二 次 名 义 
应 力 准则 。 
1.4 材料 软化 法 则 

当 笑 结 单元 应 力 状 态 达 到 损伤 初始 条 件 后 ,分 


层 损伤 将 进一步 演化 ,并 由 黏 结 单元 刚度 退化 进行 
描述 。 刚 度 退 化 由 损伤 变量 d 来 定义 ,其 计算 方法 
取决 于 不 同形 式 的 软化 准则 。ABAQUS V6.14 有 限 
元 程序 提供 了 线性 软化 法 则 与 指数 软化 法 则 ,而 两 
者 又 可 分 为 基于 有 效 位 移 和 基于 断裂 能 的 方法 进行 
运算 。 线 性 软化 法 则 与 指数 软化 法 则 的 区 别 在 于 它 
们 对 损伤 变量 d 的 计算 方式 不 同 。 在 线性 软化 法 则 
中 ,损伤 变量 d 随 着 相对 位 移 的 增加 而 线性 增长 ; 而 
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在 指数 软化 法 则 中 ,损伤 变量 d 则 随 着 相对 位 移 的 
增加 呈 指 数 增长 。 而 基于 有 效 位 移 和 基于 上 断裂 能 两 
种 软化 法 则 ,其 不 同 点 在 于 定义 损伤 变量 d 的 函数 
不 同 。 基 于 有 效 位 移 的 材料 软化 法 则 中 ,材料 的 损 
伤 变量 是 材料 损伤 发 生 后 单元 所 发 生 的 总 体位 移 的 
函数 ;基于 断裂 能 的 损伤 演化 法 则 中 损伤 变量 是 材 
料 损伤 发 生 后 材料 完全 破坏 时 所 需要 的 断裂 能 的 函 
数 。 当 损伤 变量 d 的 值 达到 1 时 ,材料 完全 破坏 。 
下 面 分 别 给 出 基于 有 效 位 移 和 基于 断裂 能 的 线性 软 
4e Yu 7 和 指数 软化 法 则 ” , 共 4 种 形式 的 表 

基于 有 效 位 移 和 基于 上 断裂 能 的 线性 软化 法 则 表 


Ela Seg EE) 

> E ACGA = 8.) e?) 
其 中: 8, 是 在 加 载 过 程 中 的 有 效 位 移 ; 8", 是 损伤 发 
生 时 的 有 效 位 移 , 由 损伤 发 生 准则 所 决定 ; 8 是 在 
对 总 模式 下 材料 完全 破坏 时 的 有 效 位 移 ,基于 有 效 
偿 稳 的 线性 软化 法 则 中 ,该 值 由 试验 得 出 ,基于 断裂 
能 的 软化 法 则 中 ,关于 o, 的 确定 方法 将 在 下 一 节 讨 


论 可 图 1 为 采用 线性 软化 法 则 的 牵引 力 - 位 移 关 系 
RAK 


Z3 
C 


图 1 线性 软化 法 则 
Fig.1 Linear softening law 


基于 有 效 位 移 的 指数 软化 法 则 表达 式 为 


0 - ài, 
J-1- 加 | 1 - exp = a 
6 d us eq eq 


] -exp(- a) 
(8) 
其 中 : a 是 与 选取 的 材料 有 关 的 无 量 纲 参数 , 它 定 义 
了 损伤 发 展 的 速率 ;exp(%) 是 指数 函数 。 
基于 断裂 能 的 指数 软化 法 则 ,得 


8., T dô 
gaj” a (9) 


ô I 
其 中 : Ta dé Fg BU 30 5 5| HZ, Ta = 
V0。 ETE +7, ; Go 为 分 层 损伤 发 生 时 的 弹性 能 ; 


Chinay iv 2E BET 
刘 志 明 ， 等 : 纤维 增强 复合 材料 混合 模式 分 层 破坏 数值 模 报 矿 光 XIV TFT sap 


图 2 为 采用 指数 软化 法 则 的 牵引 力 - 位移 关 系 示 
意图 。 


牵引 力 


s à. à. 


图 2 指数 软化 法 则 


Fig.2 Exponential softening law 


1.5 材料 破坏 准则 


在 试 件 加 载 变 形 过 程 中 , 若 损伤 变量 d = 1 , 则 
黏 结 单元 完全 破坏 ,不 再 具有 承载 能 力 。 在 单一 模 
式 傈 载 作 用 下 ,可 通过 比较 界面 的 能 量 释 放 率 与 临 
界 能 量 释 放 率 来 判断 界面 是 否 完全 分 层 。 当 d=1 
时 ,材料 完全 破坏 时 的 有 效 位 移 8 ,与 工 . 工 、 亚 型 
单一 模式 分 层 破坏 时 的 能 量 释 放 率 临界 值 Ge (i = 
1 , 开 , 亚 ) 有 如 下 关系 。 


f. 
G 0 sn 


1c — J 0 


à 
Cis | 99 65s 


od 0， , 


f 

Gy. 7 jo r,d 8, (10) 
由 上 式 可 知 ,在 单一 模式 下 , 564 为 能 量 释放 率 
Ci 的 函数 。 材 料 本 构 曲 线 下 方面 积 为 能 量 释 放 率 
C; ,在 单一 模式 荷载 作用 下 ，G; 与 材料 单一 模式 断 
裂 时 的 能 量 释放 率 G;, 相等 时 ,节点 的 有 效 位 移 即 为 
材料 在 该 模式 下 完全 破坏 时 的 有 效 位 移 , 妈 S = 
Sla o 
而 在 复合 材料 混合 模式 界面 分 层 破坏 中 ,材料 
的 完全 破坏 将 在 能 量 释 放 率 C. 达到 任 一 单一 模式 
FRIERE Ge 之 前 发 生 。 此 时 ,混合 模式 破 

坏 时 能 量 释放 率 C. 可 表达 为 
c=] 5. T «d 8,, (11) 


在 ABAQUS 有 限 元 程序 中 , 若 采 用 基于 断裂 能 
的 软化 法 则 ,需要 设置 材料 破坏 准则 判定 材料 发 生 
完全 分 层 破坏 。ABAQUS HEIE T FET ZOE] (Power 
law) Rl BK 准则 (BK law) 两 种 材料 破坏 准则 ,通过 
保证 牵引 力 -位 移 曲 线 下 的 面积 等 于 联 立 寡 指数 准 
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则 或 BK 准则 及 式 (10) 或 式 (11) 所 确定 的 能 量 释放 
KC. ,从 而 确定 材料 完全 破坏 时 的 有 效 位 移 6‘ 。 

过 指 数 准 则 的 表达 式 为 

人 GT 人 

C0 
其 中 : G1 .Gr 和 Gn 分 别 为 当前 荷载 下 I 型 本 型 和 
于 型 模式 的 能 量 释放 率 。 当 式 (12) 成 立时 , G, = 
Gr = Gi + Gy + Gno 和 窜 指 数 准 则 的 指数 yy 的 取 值 可 
以 是 1 或 2。CAMANHO 等 的 研究 表明 竹 指 数 准 
则 的 指数 y 取 值 为 1 更 适合 用 于 模拟 AS4/PEEK 2f 
维 增强 复合 材料 混合 模式 下 的 分 层 破坏 。 在 本 研究 
的 MMB 有 限 元 模型 中 ,采用 同 种 材料 进行 模拟 , 当 
材料 破坏 准则 选取 为 帘 指 数 准 则 时 , y 取 值 为 1。 

BK 准则 由 BENZEGGAGH 等 提出 ,其 表达 
AH 

G, +G\” 
= re i) | | (13) 


(12) 


J&p: c, 等 于 模型 在 当前 荷载 下 分 别 在 I 型 、 开 型 
和 加 型 分 层 破坏 中 的 能 量 释放 率 之 和 , 即 Gy = Gi + 
Gs Gy im 为 该 准则 的 材料 参数 。 文 献 [31] 给 出 
也 依据 分 层 破 坏 混合 模式 的 混合 比例 Gn/ Cr KR 
华 模 式 破坏 时 能 量 释放 率 C. 确定 参数 v 数值 的 计 


算 方法 。 对 于 ASA/PEEK 复合 材料 ,依据 试验 数据 
拟 合 到 的 参数 7 取 值 为 2. 284。 本 研究 算 例 采 用 同 
种 材料 ,因此 当 材 料 破坏 准则 选取 为 BK 准则 时 , 7 
取 值 为 2. 284。 


2 混合 模式 弯曲 模型 数值 模拟 


本 研究 建立 了 复合 材料 混合 模式 弯曲 破坏 试 件 
有 限 元 模型 模拟 复合 材料 混合 模式 分 层 破 坏 , 并 与 
试验 从 载 -位 移 曲线 进行 了 对 比 ,探讨 内 唆 力 模型 的 
材料 参数 如 界面 刚度 .单元 尺寸 .界面 法 向 和 剪 切 强 
度 .断裂 能 以 及 不 同 材 料 软化 法 则 和 材料 破坏 准则 
对 数值 模拟 准确 性 和 收敛 性 的 影响 。 

本 研究 对 文献 131] 报道 的 ASA/PEEK 复合 材 
料 层 合板 混合 模式 弯曲 破坏 试验 进行 有 限 元 分 析 ， 
试验 模型 简 图 如 图 3 所 示 。F 是 施加 的 力 ,c 是 加 载 
臂 的 长 度 ,e 的 值 与 分 层 破 坏 混 合 模式 的 混合 比例 
Gy/ G, AXR, IEAk Gr = G1 + Gr + Gn ,该 算 例 中 ， 
Gn = 0。 混 合 模 式 弯曲 试 件 长 1 = 102 mm, 宽 4 = 
25.4 mm, JE 2h =3. 12 mm, MiA REK a = 34.1 mm, 
复合 材料 单 层 与 界面 层 材料 属性 分 别 见 表 1 和 表 2。 
本 研究 取 Gr Gr = 0.5 ,此 时 ec = 44.37 mm, 


表 1 ASA/PEEK 复合 材料 层 材料 属性 
Tab.1 Material properties of AS4/PEEK composite 


"em E, /MPa E, = E, /MPa G2 = G /MPa 


G5, /MPa 


715 = Vz U53 


122 700 10 100 5500 


3 700 0.25 0.45 


表 2 界面 层 材料 参数 


Tab.2 Material properties of interface layers 


Gy. ] UN * mm^!) 


o? /MPa 


T? = 7? /MPa Km /(N - mm?) 


0.969 1.719 


80 100 1 x106 


有 限 元 模型 为 简化 建 模 省 略 了 加 载 臂 ,直接 将 


图 3 混合 模式 弯曲 试验 模型 简 图 
Fig.3 Experimental setup of the mixed-mode 


bending ( MMB) test 


MRJ F 8 4579 DESRCE TAE PARRER Fu, 和 加 
载 在 试 件 端 部 的 线 荷 载 R. ,有 限 元 模型 如 图 4 所 
示 。 而 由 于 Fa n Fe 大 小 与 加 载 辟 的 长 度 c。 有 关 , 加 
载 臂 长 度 。 又 由 破坏 模式 混合 比例 Cr Gr 所 决定 ， 
Bi FIF. 的 比例 可 由 预 设 的 破坏 模式 的 混合 比例 
直接 得 出 ， 具 体 计算 方法 可 参考 文献 [32] 。 加 载 
J F, 和 FF. 与 破坏 模式 混合 比例 的 关系 如 表 3 
所 示 。 
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| | 
z -— 
| [156mm 
EPUM T La 7] 
到 4 混合 模式 弯曲 试 件 几 何 形状 及 尺寸 图 
Fig.4 | Geometry and dimensions of MMB specimen 
本 人 研究 取 Gi” Gr = 0.5, F, = 1.87F, F, = 
0.878。 而 对 应 加 载 点 的 竖 向 位 移 w 可 由 试 件 中 点 
的 竖 向 位 移 ww 和 试 件 端点 的 竖 向 位 移 w. 通过 式 
(14) 换 算得 出 
w = Ou = 7,0, (14) 


表 3 MMB 试 件 荷 载 关系 表 
Tab.3 Mixed-mode bending ( MMB) load relationship 


Kcr 6 F F 

7 0.0( DCB) 0.00 1.00F 

Q o2 3.15F 2.15F 

e 

" 0.5 1.87F 0.87F 

cO 

Q os 1.56F 0.56; 

co 

AI 0(ENF) 1.00F 0. 00 
T = T E N O OTTA 

GE :DCB 代表 double cantilever beam 即 双 悬臂 梁 试 件 , 用 以 研究 I 
型 层 间 破坏 ;ENF 代表 end notch flexure, 即 端 部 缺口 弯曲 试 件 ,用 以 
ME RADEK. 


图 5 试 件 单元 划分 示意 图 
Fig.5 FE mesh of the MMB specimen 

在 有 限 元 模型 中 , 若 采用 的 材料 本 构 关 系 含 软 
化 段 (如 本 人 研究 的 黏 结 单元 本 构 关 系 ) ,数值 模拟 容 
易 产 生计 算 困 难 而 导致 有 限 元 分 析 提 早 终 止 。 同 时 
若 采 用 力 加 载 方式 ,将 无 法 获得 荷载 -位 移 曲 线 的 下 
降 段 ,本 研究 采用 位 移 加 载 方式 。 根 据 式 (14) 跟 简 
化 模型 (图 4) 的 几何 关系 ,将 相应 于 文献 131] 中 试 
验 曲线 终点 对 应 的 竖 疝 位 移 , 即 加 载 臂 右 端点 产生 
坚 直 向 下 6 mm 位 移 , 转 化 为 加 载 于 试 件 左 端 边线 
(图 4 所 示 .处 ) 与 中 线 位 置 处 (图 4 所 示 F.A) 的 
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位 移 ,分 别 为 向 上 6.23 mm 与 向 下 0.3 mm。 复 合 材 
料 层 采 用 SC8R 连续 壳 单 元 模拟 , 界面 层 采用 
COH3DS 黏 结 单元 模拟 。 在 下 面 所 有 数值 模拟 当 
中 ,如 图 5 所 示 , 基 于 结构 的 受 力 特性 ,复合 材料 层 
和 界面 层 沿 厚度 方向 仅 放置 一 层 单元 , 沿 宽度 方向 
放置 10 个 单元 ,在 试 件 长 度 方向 按照 材料 受 力 特 性 
分 为 3 个 分 区 ,第 一 分 区 (Region 1) WM AREK, 
此 区 域 厚度 方向 两 层 复 合 材 料 层 板 不 相连 , 层 间 不 
放置 黏 结 单元 , 且 每 层 沿 试 件 长 度 方向 划分 为 70 个 
单元 ,第 二 分 区 (Region 2) 长 度 为 试 件 长 度 一 半 与 
预 留 裂缝 长 度 之 差 , 该 区 域 在 两 复合 材料 层 板 间 放 
置 一 层 黏 结 单元 ,为 准确 预测 该 分 区 内 黏 结 单元 应 
力 分 布 , 沿 试 件 长 度 方向 划分 了 50 个 单元 ,第 三 分 
区 (Region 3) 为 试 件 右 半 部 分 ,复合 材料 层 合 板 达 
到 承载 力 峰值 直到 材料 失去 承载 能 力 , 分 层 破坏 均 
未 拓展 到 该 区 域 ,该 区 域 未 插入 黏 结 单元 并 采用 粗 
网 格 ,复合 材料 层 板 间 认 为 是 完美 蒜 结 , 沿 试 件 长 度 
方向 划分 为 100 个 单元 ,以 提高 运算 效率 。 在 下 述 
数值 模拟 中 ,如 未 特殊 说 明 , 界 面 层 材料 属性 按 表 2 
取 值 , 秋 结 单元 材料 损伤 发 生 法 则 采用 二 次 名 义 应 
力 准 则 、 材 料 软 化 法 则 采用 基于 断裂 能 的 指数 软化 
法 则 、 材 料 破坏 准则 采用 BK 准则 ,BK 准则 参数 n 
取 值 为 2. 284。 因 采用 了 含 软化 段 的 本 构 关系 ,有 
限 元 数值 模拟 容易 产生 计算 困难 而 提前 终止 ,本 研 
究 所 有 数值 计算 , 当 计 算 收 敛 效果 好 时 ,相应 于 加 载 
辟 右 端 加 载 点 加 载 位 移 达到 6 mm 时 ,计算 自动 终 
止 。 当 产生 数值 计算 困难 时 ,程序 无 法 继续 计算 而 
提前 自动 终止 或 因 计 算 时 间 过 长 而 人 为 终止 程序 继 
EAE ,计算 奏 载 -位 移 曲线 如 图 6 ~ 13 所 示 , 参 数 
化 模拟 结果 如 表 4 ~10 所 示 。 


2.1 界面 刚度 参数 化 模拟 


为 考虑 界面 刚度 值 的 大 小 对 混合 模式 弯曲 破坏 
数值 模拟 准确 性 的 影响 ,本 节选 用 一 系列 数量 级 递 
增 且 具有 代表 性 的 刚度 值 并 将 数值 模拟 结果 与 试验 
结果 进行 比较 ,为 刚度 天 | Ki 的 取 值 提供 参考 。 

所 选用 的 刚度 值 K 分 别 为 5 x 107,1 x 10*, 
6. 47 x 10* 1x10 .1.62 x10° 5 x 105,1 x 10* N/mm’, 
HB K =6. 47 x 100! N/mm? ,1. 62 x 10 N/mm 为 根 


" E y A^ bm 
据 按 公式 k- eg =20,50 计算 得 出 。 这 两 个 


取 值 可 保证 界面 变形 对 结构 整体 刚度 的 损失 小 于 
5% 和 2%。 本 小 节 共 进行 了 3 组 参数 的 模拟 :第 1 
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组 ,将 Ki 作为 变量 ,并 令 K = Kg 21 x 105N/mm!, 
预测 荷载 峰值 与 误差 ,运算 时 间 与 增 量 步 数 见 表 4。 
通过 对 MMB 模型 进行 界面 刚度 K, 参数 化 模 
拟 得 到 的 图 3 所 示 混 合 模式 弯曲 试验 模型 加 载 辟 右 
端点 的 加 载 力 (了 ) -位移 (w) 曲线 如 图 6 所 示 。 从 图 
中 可 知 , 当 界面 刚度 天 取 值 设 定 大 于 等 于 1.62 x 
10 N/ mm? 时 ,在 荷载 加 载 完毕 时 ,界面 层 并 未 发 生 
分 层 破坏 ,数值 模拟 结果 与 实际 情况 不 符合 。 当 界 
面 刚度 Ki 为 6.47 x 10* N/mm? 及 1 x 10 N/mm’ 
时 ,界面 分 层 破 坏 发 生 材料 属性 退化 ,对 应 的 丛 载 - 
位 移 曲 线 出 现下 降 段 ,其 中 ,界面 刚度 Ki1 为 1 x 
10° N/mm’? 曲线 在 达到 荷载 峰值 后 出 现 收敛 困难 ， 
运算 效率 低 ;界面 刚度 KA 6. 47 x 10* N/mm 曲线 荷 


ChinaXiv 合 作 期 刊 第 40 郑 


载 峰值 较 大 于 试验 值 ;而 界面 刚度 五 为 1 x 10* N/mm? 
曲线 下 降 段 平滑 ,预测 与 试验 曲线 接近 。 


界面 刚度 (N/mm ) 

=- - -试验 数据 

一 上 一 10 N/mm’ 
—0— 5X 10^ N/mm? 
—2— 1.62 X 10° N/mm? 
一 一 10 N/mm’? 

—€— 647 X 10* N/mm? 
—CO- 10 N/mm’ 
—fk— 5X 10! N/mm? 


S 


加 载 力 F/N 


竖 向 位 移 w/mm 


图 6 界面 刚度 Ki 参数 化 模拟 结果 


Fig.6  Interfacial stiffness K, parametric simulation results 


表 4 界面 刚度 大; 参数 化 模拟 结果 


S Tab.4 Interfacial stiffness K, parametric simulation results 

hj /(N * mm^?) 荷载 峰值 /N 与 试验 值 差 值 /% CPU 运算 时 间 /inin 增 量 步 数 备注 
试验 数据 277.0 = x Z E 
o 105 - - 13.7 100 未 分 层 
e 5 x 10? 要 = 12.7 100 未 分 层 
OO 1.62 x105 = = 165.9 397 未 分 层 
© iw 262.0 -5.42 241.9 507 收敛 困难 
x 6.47 x 10^ 322.4 16.30 85.1 220 a 
O 104 254.0 -8.30 42.4 136 - 
N 5 x 10? - z 1427 2092 不 收敛 


.全 第 2 组 , 令 Ki = Ky 21 x10* N/mm , Ky Wt 
APERIRA EXER MMB 模型 进行 界面 刚度 K 参 


v GS Ha php AN ` — 
数 屁 模拟 得 到 的 力 -位 移 曲 线 如 图 7 所 示 。 

ge UA ge mem 3 

界面 刚度 (N/mm) 
= -- -试验 数据 
© 一 一 10 Nmm 

—0—5X 10! N/mm' 
—^— 1.62 X 10° N/mm* 
一 一 10 N/mm* 
—0— 5X 10 Nmm 
一 一 10 N/mm? 
—*— 5X 10° N/mm 


加 载 力 F/N 


0 1 2 3 4 5 6 
竖 向 位 移 w/mm 


K7 界面 刚度 Ki 参数 化 模拟 结果 
Fig.7 Interfacial stiffness K, parametric simulation results 
从 图 中 可 知 , 随 着 界面 刚度 Ki 的 降低 ,曲线 斜率 
随 之 降低 。 且 当 界 面 刚 度 K1 取 值 1 10* N/mm ,只 
改变 界面 刚度 Ki 的 取 值 时 ,由 于 界面 刚度 K| 取 值 
过 大 ,在 荷载 加 载 完全 时 未 发 生 界 面 分 层 , Ky 的 取 
值 大 小 只 对 曲线 斜率 有 影响 ,曲线 并 未 出 现下 降 段 ， 


此 组 分 析 未 出 现 分 层 破坏 。 
由 于 复合 材料 界面 混合 模式 分 层 破坏 包含 型 
及 开 型 分 层 破坏 , 故 下 面 同时 令 Ki、Ki 为 变量 , 取 前 
述 一 系列 刚度 值 ,并 取 Ki = Ky; Kg 21 x 10* N/mm’, 
预测 荷载 峰值 与 误差 ,运算 时 间 与 增 量 步 数 见 表 5。 
第 3 组 ,通过 对 MMB ILT K, Ky 刚度 参 
数 化 模拟 得 到 的 力 -位 移 曲线 如 图 8 所 示 。 
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界面 刚度 (N/mm) 
-一 试验 数据 
一 一 10 N/mm’ 
—0— 5X 10 N/mm’ 
—£8— 1.62 X 10 N/mm? 
一 一 10 N/mm? 

—9— 6.47 X 10° N/mm? 
一 0 一 10 N/mm? 
—&— 5X 10! N/mm? 


竖 向 位 移 w/mm 


FR] 


8 ERBEK, Kı 参数 化 模拟 结果 
Fig.8 Interfacial stiffness K; ,and Ky 


parametric simulation results 


从 图 中 可 知 当 刚度 取 值 大 于 等 于 5 x 10 N/m’ 
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时 ,界面 层 并 未 发 生 分 层 破 坏 ,数值 模拟 结果 与 实际 
情况 不 符合 。 当 界面 刚度 K， = Kg 分 别 为 5 x10 、 
1 x10* ,6.47 x10* ,1 x 10^ 1.62 x10 N/mm’ 时 ， 
界面 分 层 破 坏 发 生 材 料 属性 退化 ,对 应 的 荷载 -位 移 
曲线 出 现下 降 段 ,其 中 ,界面 刚度 为 1. 62 x 10 N/mm? 
及 1 x10 N/ mm? 曲线 在 荷载 达到 峰值 后 小 幅 又 
降 ; 而 界面 刚度 为 6.47 x 10* N/mm ,曲线 在 荷载 达 
到 299. 5 N 时 ,由 于 分 层 损 伤 发 展 迅 速 , 黏 结 单元 刚度 
退化 ,结构 承载 力 下 降 ,曲线 出 现下 降 段 , 答 载 峰值 较为 
接近 试验 值 ;界面 刚度 分 别 为 1 x 10 "及 5 x10” N/mm’ 
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的 曲线 在 荷载 达到 峰值 后 均 发 生 震 荡 , 且 并 未 出 现 
下 降 段 ,与 实际 情况 不 符合 。 

上 述 分 析 绪 果 表 明 ,界面 刚度 K, 主要 影响 曲 
线 峰值 ,而 界面 刚度 Kr 主要 影响 曲线 斜率 ;数值 模 
拟 的 力 - 位 移 曲 线 的 形状 ,由 界面 刚度 K Ka 共同 
影响 。 当 界面 刚度 K,. Kg 共同 取 值 按 公 式 K = 


Bs ,并 取 有 为 20 计算 得 到 的 值 6.47 x 10° N/mm 
时 数值 模拟 结果 较为 准确 , 且 曲 线 较为 平滑 ,下 述 研 
究 将 采取 此 刚度 值 进行 数值 模拟 。 


RS RANEK, Kı 参数 化 模拟 结果 


Tab.5 Interfacial stiffness K, and Ky parametric simulation results 
Ki Kq/](N * mn) 荷载 峰值 /N 与 试验 值 差 值 /% ”CPU 运算 时 间 /min 增 量 步 数 备注 
试验 数据 277.0 i z - i 
10$ - - 8.1 100 未 分 层 
im 5 x 10? - - 3.2 102 未 分 层 
T= 1.62 x 10 324.5 17.1 335.2 680 误差 较 大 
335.7 21.2 93.3 234 误差 较 大 
299.5 8.1 38.5 240 较 准 确 
247.1 -10.8 279.0 533 RER 
AT 5x10 258.9 -6.5 311.9 610 发 生 震 荡 


20! 单元 尺寸 参数 化 模拟 
:二 因为 内 聚 区 里 应 力 变化 梯度 大 ,为 得 到 准确 的 


别 勇 0.1.0.25 .0.338( Balzani 9X fi P" ,将 长 为 
16:9 mm 的 范围 内 划分 为 50 个 单元 ) ,0.5、1 mm, 
图 9 所 示 为 不 同 单元 尺寸 的 复合 材料 混合 模式 分 层 
破坏 数值 模拟 得 到 的 荷载 -位 移 曲 线 ,相关 模拟 数据 
见 表 6。 


单元 尺寸 /mm 
--- 试验 数据 
== | mm 

— — 0.5 mm 
—— 0.338 mm 


e 
N 
w 
AR 
[i 
CN 


竖 向 位 移 w/mm 


图 9 单元 尺寸 参数 化 模拟 结 


Fig.9 Mesh size parametric simulation results 


表 6 各 单元 尺寸 模型 预测 荷载 峰值 与 CPU 运算 时 间 
Tab.6 Predicted peak loads and computational time 


of FE models with different mesh sizes 


单元 尺寸 /mm 荷载 峰值 /N CPU 运算 时 间 /min 
试验 数据 277.0 = 
1 无 下 降 段 18.5 
0.5 285 23.9( 收敛 困难 ) 
0.338 299.5 106.5 
0.25 298.2 143.1 
0.1 277.8 371.3 


结果 显示 ,采用 较 大 的 单元 尺寸 (1 mm ) 得 到 的 
模拟 曲线 并 未 出 现下 降 段 ,与 试验 曲线 并 不 吻合 :其 
余 采 用 不 同 单元 尺寸 的 数值 模拟 得 到 的 预测 曲线 ， 
其 荷载 峰值 大 小 根据 选取 的 单元 尺寸 不 同 而 有 所 不 
同 。 当 单元 尺寸 采用 0.5 mm 时 ,曲线 达到 荷载 峰值 
(285.0 N) 后 出 现 收敛 困难 , 而 随 着 单元 尺寸 的 减 
小 ,数值 模拟 荷载 峰值 逐渐 接近 试验 值 。 当 单元 尺 
才 采 用 0. 25 mm 时 ,数值 模拟 荷载 峰值 为 298.2 N, 
与 单元 尺寸 采用 0.338 mm 时 和 荷载 峰值 (299.5 N) JH 
比 , 略 接近 试验 值 。 当 单元 尺寸 采用 0. 1 mm 时 ,和 荷 
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载 峰值 为 277. 8 N, 允 近 和 荷载 峰值 。 从 运算 效率 的 角 
度 考虑 , 随 着 单元 尺寸 的 减 小 ,单元 数量 较 多 ,计算 
效率 较 低 , 当 单 元 尺寸 采用 0. 1 mm 时 ,CPU 运算 时 
间 为 371.3 min, 所 耗费 时 间 为 单元 尺寸 采用 
0.338 mm 时 的 约 3.5 fi. 

综 上 ,采用 较 小 单元 尺寸 (单元 尺寸 小 于 
0.5 mm) 的 黏 结 单元 模拟 复合 材料 混合 模式 分 层 破 
坏 准 确 性 较 高 , 当 单 元 尺寸 过 小 时 (单元 尺寸 为 
0.1 mm) ,由 于 单元 数量 相对 增多 ,计算 效率 低 。 
此 从 计算 效率 和 精度 的 角度 综合 考虑 ,在 后 述 的 研 
究 中 ,单元 尺寸 取 0.338 mm, 


2.3 界面 强度 参数 化 模拟 


一 本 节 探 讨 界 面 层 强度 对 数值 模拟 的 影响 , 进 

了 油 组 参数 的 模拟 :第 1 组 ， r o° 
HD0, 40,60,80,100 MPa, 剪 切 强度 取 100 MPa( 参 
ZIE) ;第 2 组 , 令 界面 层 剪 切 强度 W 20,40,60, 
80100 MPa ,法 向 强度 取 80 MPa( 参考 值 ) 。 在 长 度 
并 条 上 秋 结 单元 尺寸 为 0. 338 mm。 数 值 模拟 结 
Aff] 10 和 图 11 所 示 , 不 同 界面 法 向 强度 及 剪 切 强 
庆 对 应 的 预测 荷载 峰值 见 表 7, 表 8。 


界面 法 向 强度 /MPa 
-一 -试验 数据 
—L1- 100 MPa 
—QO- 80 MPa A 
—3jdÀ— 60 MPa x 
--/x- 40 MPa 


竖 向 位 移 w/mm 


图 法 向 强度 参数 化 模拟 结果 


Interfacial normal strength parametric 


图 10 界 
Fig. 10 


simulation results 


界面 剪 切 强度 /MPa 
300} - - -试验 数据 
一 100 MPa 
250} —O— 80 MPa 
zZ 一 大- 60 MPa 
E 200F --&- 40 MPa 
ps —V— 20 MPa 
& 150 
100 
50 
0 1 2 3 4 3 6 
竖 向 位 移 wymm 
图 11 界面 剪 切 强度 参数 化 模拟 结果 
Fig. 11 Interfacial shear strength parametric simulation results 
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表 7 界面 法 向 强度 参数 化 模拟 结果 


Tab.7 Interfacial normal strength parametric simulation results 


on /MPa ”荷载 峰值 /N 与 试验 值 差 值 /% CPU 运算 时 间 /min 
试验 数据 277.0 = z 

100 302. 1 9.1 92.5 
80( 参 考 值 ) 299.5 8.1 106.5 

60 302.9 9.4 70.3 

40 300.5 8.5 63.1 

20 288.2 4.0 129.1 


表 8 界面 剪 切 强度 参数 化 模拟 结果 


Tab.8 Interfacial shear strength parametric simulation results 
TP/MPa FIREAN 与 试验 值 差 值 /% CPU 运算 时 间 /min 
试验 数据 277.0 m 

100( 参 考 值 ) 299.5 8.1 106.5 
80 300.7 8.6 87.6 
60 299.9 8.3 87.6 
40 299.9 8.3 75.6 
20 295.9 6.8 138.0 


数值 模拟 结果 表明 ,界面 层 极限 强度 的 取 值 对 
数值 模拟 结果 影响 很 小 。 随 着 界面 强度 的 下 降 ,各 
预测 荷载 峰值 十 分 接近 ,荷载 峰值 预测 值 与 试验 数 
据 大 致 接近 。 另 外 , 随 着 界面 强度 的 下 降 , 除 界面 强 
度 选取 20 MPa 的 数值 模拟 外 ,CPU 运算 时 间 也 随 之 


减 小 。 对 于 混合 模式 分 层 破坏 的 模拟 结果 而 言 , 界 
面 强度 适当 地 降低 带 来 的 影响 为 荷载 峰值 的 略微 降 


低 ,界面 强度 的 提高 则 会 引起 数值 模拟 运算 时 间 的 
增加 。 综 上 ， e MATES RE 
面 法 向 强度 ,提高 运算 效率 。 


2.4 断裂 能 参数 化 模拟 


本 节 分 析 中 ,采用 基于 断裂 能 的 指数 软化 法 则 ， 
震 指 数 准 则 作为 破坏 准则 。 根 据 式 (12) 可 知 ,在 复 
合 材 料 发 生 混合 模式 分 层 破坏 时 ,其 断裂 能 G 的 大 
小 ,与 所 选取 的 I 型 破坏 断裂 能 Ci. 及 工 型 破坏 断 
ABE GI, 相关 。 因 此 ,对 黏 结 单元 而 言 ,其 本 构 关 系 
曲线 下 面积 即 等 于 断裂 能 C. 的 大 小 , 即 本 构 关 系 曲 
线 受 Ci Gr: 取 值 所 影响 。 为 探讨 断裂 能 对 MMB 
模型 数值 模拟 结果 的 影响 ,选用 界面 刚度 K, = 天 
6.47 x 10 N/mm , Ky = 1 x 10* N/mm 界面 强度 
o? = 80 MPa, T° = 100 MPa 界面 层 断 裂 能 Ci 及 
Gie 同时 设置 为 实际 值 的 0.8、0.9、1、1.1、1.2 fis 
其 荷载 -位 移 曲线 如 图 12 所 示 ,和 荷载 峰值 及 模型 运 
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算 时 间 如 表 9 所 示 。 


断裂 能 
- -试验 数据 
0.8 倍 


竖 向 位 移 w/mm 
图 12 断裂 能 参数 化 模拟 结果 
Fig. 12 Fracture energy parametric simulation results 
ROI MMB 模型 断裂 能 参数 化 模拟 结果 
Tab.9 MMB model fracture energy parametric 


y= simulation results 


BARA ”荷载 峰值 /N SRRA CPU 运算 时 间 /min 


m——— 271.0 = - 


ces - = 112.5( 不 收敛 ) 
C.» 无 下 降 段 - 169.5 
) 299.5 8.1 38.5 
311.9 9.3 75.86 
无 下 降 段 一 84.7 


AM 数值 模拟 结果 表明 ,对 复合 材料 混合 模式 分 层 


T E ,断裂 能 的 改变 对 荷载 峰值 及 对 应 数值 模 
拟 申 线 的 形状 影响 十 分 大 。 若 断裂 能 取 值 与 实际 值 
相差 过 大 ,不 仅 会 影响 结构 的 承载 力 ,而 且 对 应 的 预 
测 圭 线 将 与 实际 不 符 。 作 为 试验 测 得 的 材料 参数 ， 
界面 层 断 裂 能 的 取 值 必须 足够 准确 以 保证 数值 模拟 
的 锥 确 性 。 


2.5 材料 软化 法 则 及 材料 破坏 准则 参数 化 模拟 


ABAQUS V6. 14 有 限 元 程序 对 黏 结 单元 在 材料 
软化 阶段 的 材料 行为 提供 了 基于 位 移 和 基于 断裂 能 
的 线性 和 指数 软化 法 则 。 材 料 软化 法 则 描述 了 材料 
在 发 生 损伤 后 ,牵引力 - 位 移 曲 线 的 下 降 方式 。 选 取 
不 同 的 材料 软化 法 则 ,将 影响 竹 结 单元 的 本 构 关 系 
曲线 的 形状 ,对 复合 材料 在 分 层 破坏 过 程 中 力学 行 
为 的 定义 具有 重要 意义 。 本 节 将 讨论 选用 基于 有 效 
位 移 和 基于 断裂 能 的 线性 和 指数 软化 法 则 对 运算 结 
果 的 影响 。 

下 面 将 对 比 基 于 有 效 位 移 和 基于 断裂 能 的 线性 
软化 法 则 和 指数 软化 法 则 (其 中 ,基于 断裂 能 的 软化 
法 则 考虑 寡 指 数 准 则 和 BK 准则 两 种 不 同 材 料 破坏 准 


Chinay nr BB3T 
刘 志 明 ， 等: AERE A EC A Bo e cott IV TERR TI sss 


则 的 影响 ) , 共 6 组 数值 模拟 在 运算 结果 上 的 差异 。 

本 研究 取 Cr Cr = 0.5 ,相应 地 材料 断裂 能 G. 
取 1.131 N/mm ^ ,对 于 基于 有 效 位 移 的 线性 软化 
法 则 ,材料 完全 破坏 时 的 有 效 位 移 8&, =2 CT HC 
Te - (05). * (r5). * (00). ,得 8 20.0139 mm; 对 
于 基于 有 效 位 移 的 指数 软化 法 则 ,通过 联 立 公式 
(1)、(2)、(8) 和 (10) 计 算得 出 G1AGr =0 时 , 工 型 
破坏 位 移 8 = 0.00455 mm, Gy/ Gr = 1 时 ,本 型 
破坏 位 移 8 = 0.007 77 mm, 有效 位 移 8 由 
ABAQUS 根据 51 及 61 确定 。 数 值 模拟 结果 如 图 13 
所 示 , 相 关 计算 数据 如 表 10 所 示 。 

由 图 13 和 表 10 可 看 出 ， 基 于 位 移 的 软化 法 则 
与 基于 断裂 能 的 软化 法 则 在 数值 模拟 结果 上 差别 很 
大 。 基 于 位 移 的 指数 软化 法 则 在 位 移 加 载 到 
2.6 mm 时 达到 荷载 峰值 135. 5 N, 在 位 移 加 载 到 
4.1 mm 计算 模型 出 现 收敛 困难 的 情况 ; 而 基于 位 移 
的 线性 软化 法 则 在 位 移 加 载 到 3.5 mm 时 达到 荷载 
峰值 195.7 N ,与 试验 值 相差 较 大 。 采 用 基于 断裂 能 
的 线性 软化 法 则 和 和 窜 指 数 准则 作为 材料 破坏 准则 的 
分 析 , 荷 载 在 达到 262.0 N 后 曲线 出 现 轻微 震荡 ,和 茶 
载 峰值 为 269.5 N ,与 试验 值 较 接 近 ; 基 于 断裂 能 的 线 
性 软化 法 则 和 BK 准则 作为 材料 破坏 准则 的 分 析 , 预 
测 曲 线 在 达到 荷载 峰值 242. 1 N 后 亦 出 现 轻微 震荡 。 

ibis e E RAN 


IBrA Re B Ze PEKALA 0 — n6 48 AE 
fek 多 化 法 则 一 一 BK 准则 


300 [一 一 基于 位 移 的 线性 软化 法 则 
:于 位 移 的 指数 软化 法 则 


竖 向 位 移 w/mm 
图 13 材料 软化 法 则 及 材料 破坏 准则 参数 化 模拟 结果 


Fig.13 Material softening law and delamination propagation 


criterion parametric simulation results 

采用 基于 断裂 能 的 指数 软化 法 则 和 BK 准则 作 
为 材料 破坏 准则 的 分 析 , 在 位 移 加 载 到 5.06 mm 时 
得 到 的 预测 曲线 荷载 峰值 (299. 5 N) 与 试验 值 
(277.0N) 较 为 接近 。 采 用 基于 断裂 能 的 指数 软化 
法 则 和 有 需 指 数 准 则 作为 材料 破坏 准则 的 分 析 在 位 移 
加 载 到 5. 31 mm 时 曲线 达到 了 和 荷载 峰值 326.3 N , 较 
大 于 试验 值 。 在 运算 效率 方面 ,采用 基于 断裂 能 的 
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指数 软化 法 则 和 短 指 数 准 则 作为 材料 破坏 准则 的 分 
析 所 需要 的 CPU 运算 时 间 为 210.3 min ,而 采用 基于 
断裂 能 的 指数 软化 法 则 和 BK 准则 作为 材料 破坏 准 
则 的 分 析 的 运算 时 间 为 106.5 min。 可 以 看 出 ,选用 
基于 断裂 能 的 指数 软化 法 则 和 BK 准则 作为 材料 破 
坏 准 则 能 更 好 地 提高 运算 精度 和 运算 效率 。 


数值 模拟 表明 ,采用 不 同 的 软化 法 则 及 材料 破 
坏 准则 影响 数值 模拟 的 收敛 性 和 效率 。 选 用 基于 断 
裂 能 的 指数 软化 法 则 能 准确 地 模拟 复合 材料 混合 模 
式 分 层 破坏 ,其 中 选用 BK 准则 作为 材料 破坏 准则 
的 计算 收敛 性 展 好 ,数值 计算 取得 较 好 的 结果 , 且 计 


表 10 软化 法 则 及 材料 破坏 准则 参数 化 模拟 结果 


Tab. 10 Material softening law and delamination propagation criterion parametric simulation results 


软化 法 则 荷载 峰值 人 N 与 试验 值 差 值 /% CPU 运算 时 间 /min 增 量 步 数 
试验 数据 277.0 - = = 
基于 断裂 能 的 线性 软化 法 则 一 一 晖 指数 准 则 269.5 -2.7 80.3 370 
基于 断裂 能 的 线性 软化 法 则 一 一 BK 准则 242.1 -12.6 136.0 710 
AE T Br fe KITE C Ls M—F 8 28 EE DUI 326.3 17.8 210.3 815 
基于 断裂 能 的 指数 软化 法 则 一 一 BK 准则 299.5 8.1 106.5 240 
基于 位 移 的 线性 软化 法 则 195.7 -29.4 254.1 961 
基于 位 移 的 指数 软化 法 则 135.5 -51.1 413.4 2 086 
~ 5) 合 理 地 选用 正确 的 材料 软化 法 则 及 材料 破坏 
3 三 结 it 准则 对 运算 效率 提高 尤为 关键 。 选 用 基于 断裂 能 的 


旺 本 研究 探讨 内 聚 力 模型 的 材料 参数 如 界面 刚 


旨 数 软化 法 则 和 采用 BK 准则 作为 材料 破坏 准则 的 
分 析 运 算 效 率 高 ,运算 结果 收敛 且 准 确 ,适用 于 模拟 


谋 ( 界 面 法 向 和 剪 切 强度 、 断 裂 能 以 及 单元 尺寸 .不 
同 烤 料 软化 法 则 和 材料 破坏 准则 等 对 复合 材料 层 合 
板 混合 模式 分 层 破坏 弯曲 试验 的 数值 模拟 准确 性 和 
收敛 性 的 影响 。 

1) 在 界面 刚度 参数 化 模拟 中 ,采用 了 不 同 界面 
刚 庆 值 进行 了 数值 模拟 ,所 需要 的 运算 时 间 以 及 有 
限 元 分 析 所 进行 的 增 量 步 数 有 较 大 区 别 ,合理 选用 
界面 单元 刚度 值 能 大 大 减少 运算 时 间 , 提 高 运算 效 


层 合板 混合 模式 分 层 破坏 。 老 指数 准则 和 BK 准则 
作为 判定 忒 结 单元 是 否 完全 破坏 的 重要 准则 ,可 按 
具体 情况 选用 适当 的 材料 破坏 准则 。 

本 研究 结果 为 更 好 地 在 有 限 元 分 析 中 采用 内 聚 
力 模型 和 务 结 单元 模拟 实际 工程 中 复合 材料 层 合板 
及 其 结构 更 为 普遍 的 分 层 破 坏 模式 , 即 混合 模式 分 
层 破 坏 的 力学 行为 提供 了 参考 指导 作用 ,有 利于 复 
合 材料 研究 人 员 或 结构 设计 工程 师 更 好 地 选择 模型 


党首 过 分 析 结果 可 知 ,对 于 复合 材料 混合 模式 界面 
AVE EK GU TURON 建议 公式 K = P^ , 取 p = 


20 时 的 刚度 值 为 界面 刚度 K: Kr 值 ,其 数值 模拟 
结果 较为 准确 。 

2) 界 面 法 向 和 剪 切 强度 的 参数 化 分 析 表 明 其 对 
混合 模式 分 层 破坏 数值 模拟 结果 影响 很 小 , 当 数 值 
计算 出 现 困 难 时 ,可 适当 降低 界面 强度 ,提高 运算 

3) 断裂 能 作为 材料 本 丑 的 性 质 之 一 直接 与 材料 
的 承载 力 相关 ,作为 试验 测 得 量 其 值 必须 足够 准确 
保证 数值 模拟 的 准确 性 。 

4) 单 元 尺寸 的 选择 影响 计算 结果 的 准确 性 和 计 
算 效 率 , 其 值 必须 足够 小 以 保证 计算 精度 和 效率 。 
当 单 元 尺寸 过 大 时 ,计算 结果 不 准确 ;过 小 时 ,计算 


参数 与 准则 以 保证 数值 计算 的 准确 性 、 收 敛 性 和 计 
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